Hysteresis in driven disordered systems: from plastic depinning to magnets 
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We study the dynamics of a viscoelastic medium driven through quenched disorder by expanding 
about mean field theory in 6 — e dimensions. The model exhibits a critical point separating a 
region where the dynamics is hysteretic, with a macroscopic jump between strongly pinned and 
weakly pinned states, from a region where the sliding state is unique and no jump occurs. The 
disappearance of the jump at the critical point is described by universal exponents. As suggested 
in [|l], the model appears to be in the same universality class as the zero-temperature random field 
Ising model of hysteresis in magnets. 



I. INTRODUCTION 



Since the eighties the depinning of various elastic inter- 
faces, ranging from domain walls in magnets to geological 
faults, has received great attention. Many elastic inter- 
face models with dissipative dynamics display a continu- 
ous transition from pinned to sliding at a critical driving 
force Fc |@,||. Near this transition the behavior is uni- 
versal on long length scales, i.e. it does not depend on 
microscopic details, but only on general properties of the 
system, such as symmetries, dimensions, range of inter- 
actions, and the dynamics. The sliding state is unique 
and no hysteresis can occur |9| . This class of systems has 
been studied extensively by renormalization group (RG) 
methods |^ and numerical simulations Meanwhile, 
a large body of experimental |^,^ and numerical work 
iQ,^ has shown that many extended condensed matter 
systems with strong disorder exhibit a spatially inhomo- 
geneous plastic response, when set into motion by an ex- 
ternal drive In this plastic flow regime, topological 
defects proliferate and the system is broken up in fluid- 
like regions flowing around pinned solid regions. The 
elastic restoring forces are replaced by viscous flow on 
various scales allowing for hysteretic response. Ex- 
amples include vortex arrays in disordered superconduc- 
tors, charge density waves (CDWs) in anisotropic met- 
als, colloids, and many others. Only for weak disorder 
these extended systems can be described as elastic ob- 
jects pulled through a quenched random medium by a 
uniform force F . Although several mean-field (MF) mod- 
els of dissipative dynamics with locally underdamped re- 
laxation have been proposed in the literature in various 
physical contexts and shown to exhibit hystere- 

sis, very little is known analytically about the behavior 
of such models in finite dimensions. This article provides 
first analytical (RG) results for a viscoelastic model that 
was introduced and previously studied in MF theory by 
Marchetti and collaborators 0|, and makes an impor- 
tant contribution towards the identification of universal- 



ity classes for nonequilibrium disordered systems. In |^ 
it was shown that in MF theory the model has a criti- 
cal point, separating continuous from hysteretic dynam- 
ics. In this paper we show that in finite dimensions the 
behavior maps onto critical hysteresis in disordered mag- 
nets as modeled by the zero temperature non-equilibrium 
random field Ising model (RFIM) |Q. By showing that 
the two models are in fact in the same universality class 
we obtain the critical exponents for hysteresis in plastic 
flow in an expansion in 6 — e dimensions. 



II. THE DRIVEN VISCOELASTIC MODEL 

In Ref. it is proposed that a description of shear de- 
formations in the plastic regime may be obtained by fo- 
cusing on the dynamics of coarse-grained degrees of free- 
dom (the solid-like regions) that are allowed to slip past 
each other. The model of a driven viscoelastic medium 
is then obtained by replacing the elastic couplings of 
displacements in the coarse-grained model of an over- 
damped driven elastic medium with Maxwell-type cou- 
plings of velocities Considering, for simplicity, the 
overdamped dynamics of a scalar field, the equation of 
motion for the local displacement at discrete lattice sites 
«, u(r, t) Ui{t), is 



= V / ds^i,je'^'-'^'^[uj{s)-u,{s)\+F + F,{ui), 

(1) 



where the dot denotes a time derivative, the summation is 
restricted to nearest neighbor pairs and 70 is the friction. 
If all the nearest-neighbor elastic couplings are equal 
(/iy = /i > 0), the first term on the right hand side of Eq. 
(|l|) is the discrete Laplacian in d dimensions. The second 
term is the homogeneous driving force, F , and Fi(ui) de- 
notes the pinning force arising from a quenched random 
potential, Vi{ui), Fi{ui) = -dVi/dui = h,Y{u, - (3i), 
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with Y{u) a periodic function with period 1 and f3i ran- 
dom phases uniformly distributed in [0,1]. The ran- 
dom pinning forces hi are distributed with probabihty 
p{h) = e~'^/'^°/ho, of width ho. The model can be used 
to describe pure shear deformations, where the interac- 
tions among the degrees of freedom are transverse to the 
direction of mean motion. As shown in the integro- 
differential equation can be transformed into a second 
order differential equation, given by 

j-T—^jii., = -u.j] + F + F.i{ui), (2) 

iw 

with rjij = fXijT. The MF critical line 77c (t) separating 
continuous from hysteretic dynamics was obtained in 0] 
and is shown in the inset of Fig. |l|. 
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where (...) denotes an average over disorder and T{hi) is 
the period of the i-th degree of freedom (i.e. the time 
over which the displacement Ui is incremented by one). 




for h<2g , 



T{h) = 00 , for h>2g , (5) 

with g = F+rfv. The various types of dynamical response 
are shown in Fig. 1. For strong coupling (relative to dis- 
order), rj > rjc, a marginally pinned degree of freedom 
pushes over its neighbors, driving a vertical hysteretic 
jump in the iJ — F curve between strongly pinned (slowly 
moving) and weakly pinned (fast moving) states. For 
weak coupling (j] < rjc) most degrees of freedom advance 
independently and the v — F curve has no macroscopic 
jumps. In MF the hysteresis disappears at (and below) a 
special value r] = ijc, where the v — F curve has a vertical 
slope. Preliminary results indicate the same for simu- 
lations in three dimensions (up to fluctuations). Figure 
^ shows a schematic phase diagram for the dynamical 
model defined by Eq. (|) in the (F, r]) plane. The MF 
critical point {F^rjc) is obtained from x = dv/dF 00 
and d^v/dF"^ 00, where x is the static response func- 
tion. As shown in Appendix A, near this critical point 
the velocity v — Vc, with Vc ~ v(r]c, Fc), obeys a universal 
scaling law, 

^--c-M'^e±(r^), (6) 



FIG. 1. Typical MF v - F curves for r = and rj = 24 (dashed-dotted), rj ^ rjc = 9.745 (solid), and r? = 6 (dashed). The 
inset shows the phase diagram in the (r/, r) plane. The shaded amahigrfehie m^n-^\i^^^^^giiypSLmi^^iiff^^^i^i..-\- rjc) , with 

/ = F — Fc{ric), is a (non-universal) rotation between the 



In this paper we consider the limit t = 0, with rj = /it ^ 
0, referred to as the viscous limit, and study the critical 
point along the t — line. The MF analysis suggests 
that that any finite value of r may be irrelevant and the 
large scale behavior be generically described by r = 0. 
In this limit the model reduces to 



+ F + F,{u,). 



(3) 



We have introduced dimensionless variables by incorpo- 
rating 7o in our unit of time and scaling all forces and 
velocities by the width ho of the disorder distribution. 
The relevant dimensionless tuning parameters are then 
the strength i] of the viscous coupling and the driving 
force, F. The MF solution of Eq. (^) for piece-wise 
parabolic pinning potential was obtained in and is 
given by 



\Tih,) 



(4) 



control variables (r, /) and the scaling variables {r,f'). 
The ± refers to the sign of r. In MF theory the critical 
exponents have the values Pmf — 1/2 and Smf — 3. 
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FIG. 2. Schematic MF phase diagram and RG flows, (a) The bold ci 
vertical jump in the MF response (i.e. where dv/dF 00) when the forc€ 
configuration at F = (up); and (ii) decreased adiabatically from -l-oo (d 
these lines, (r^c, fc(??c)). (b) RG flows below 6 dimensions in the (x~^,u) 
d = 6 and the critical behavior is obtained by linearizing around the Wils 
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III. BEYOND MEAN-FIELD THEORY 

A. Formalism 

Here we set up the formalism to study fluctuations 
about MF theory. We use the method of Martin, Siggia 
and Rose (MSR) to transform the stochastic equation of 
motion (||) into a field theory. This is done by introduc- 
ing a generating functional Z, which contains the sum 
over the probabilities of all possible paths {i.e. micro- 
scopic configurations) which the system follows for dif- 
ferent configurations of disorder as the force F is slowly 
increased. The generating functional is written as the ex- 
ponential of an action, which is then renormalized pertur- 
batively by standard methods. The MSR formalism was 
first adapted to driven systems by Narayan and Fisher to 
describe the depinning of an elastic model of CDWs [D . 
It was later employed by Dahmen and Sethna to study 
the hysteretic response of disordered magnets . Here 
we only give a very sketchy summary of the method and 
refer to the literature for details. The crucial point is 
to recognize that the coarse grained local velocities {vi 
defined below) in the present model play the role of local 
magnetizations in the RFIM. The formal manipulations 
then follow closely those of Ref. |^^, although impor- 
tant differences occur in the detailed calculation of the 
response and correlation functions. 



The functional Z is defined as a product of (5-functions 
(one for each site), each of which imposes that the dy- 
namics of each Ui is described at all times by the equation 
of motion (^. The 5- functions are then rewritten using 
a standard identity, so that 



1 = Z = [du][du]J[u]e' 



(7) 



with Ui MSR auxiliary fields. Here J[u] is a functional 
Jacobian chosen so that Z integrates to unity and S is 
an action given by 



E^»j' 



[u-j ~ ill] ~ F - Fi{ui) 



(8) 



To study perturbatively fluctuations about MF theory we 
change variables from the Ui and Ui to the local veloci- 
ties Vi = (1/r?) '7u'"i by introducing another auxiliary 
Vi. The MF solution (uo,wo) is the saddle point of 



field, 



the resulting generating functional. By expanding about 
this solution and shifting the definition oi v io v — vq^ 
(and absorbing a factor i in the definition), we obtain 
the generating functional, 



[dv] [dv\e 



•Sett 



(9) 



with an effective action 



I (^1 •)••••, ^m+n )Vj{ti)...Vj{tm)Vj{t„i+i)...Vj{tjn+n)- (10) 



The Umn are the linear and nonlinear local response and 
connected ()c) correlation functions, of the form 

d d 



B. Results 



-{u{ti)...u{t,n))u, 



(11) 



de{tm+i) de{tjn+n) 

The suffix I in Eq. ( pT| ) indicates that the response and 
correlation functions are obtained by solving the MF 
equation with fixed local fields vq and an infinitesimal 
perturbation, rie(t,ti), at ti < t, 

{l+q)dtuj{t) = 



The effective action (HO 



7m.j+F.j{uj)+r]e{t,ti). (12) 

has the same structure as that of 
the spin model in Ref. [ 14[ , with local velocities Vi replac- 
ing local magnetizations. Note that while in the RFIM at 
fixed field all spins are fixed in time (in the steady state), 
in the CDW model each phase is subject to a quasiperi- 
odic noise due to its coupling to the neighboring sites 
that move in a periodic potential even at fixed force. In 
the calculation done here this additional noise is aver- 
aged out and does not seem to affect the scaling results 
on long length scales. A detailed study whether or not 
it matters on long length scales in ways not captured by 
this approach is left to future work. 



The response function of order n is calculated by per- 
turbing the system with n (5-function pulses of strength 
£1,62, ...,£« at ti < t2 < ... < tn, and evaluating the re- 
sponse at a time t > tn. In contrast to the case of driven 
CDWs I), terms of order - e-^'^C'', with A = h/{l + ri), 
cannot be neglected in the calculation, as, even when the 
mean velocity is very small, a fraction of degrees of free- 
dom may be sliding freely, with corresponding large pe- 
riods. Furthermore, the response must be evaluated over 
the entire history to obtain the correct low frequency be- 
havior. An example of such a calculation for the simplest 
case of the linear response function is shown in Appendix 
B. Although the details of the calculation for the present 
model are different from the spin model, we find that 
the Umn of the two models have indeed the same long 
wavelength and low frequency behavior, i.e. the effective 
actions of the two models are identical. Consequently, 
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the RG analysis of the viscoelastic model is the same as 
for the RFIM iQ and is not reproduced here. (We obtain 
the same e-expansion to all orders in e. Note that while 
the e-expansion for the equilibrium RFIM is controversial 
[ prt , and the e-expansion for the nonequilibrium RFIM 
has been mapped onto that of the equilibrium RFIM 
to all orders in e, our main result here does not depend 
on the outcome of this controversy: the mapping of the 
plastic CDW depinning model onto the hysteretic RFIM 
is based on mapping the actions of the two models onto 
each other, by identifying the local velocities in the CDW 
model with the local magnetizations in the spin model. 
This mapping holds independently of the subsequent e- 
expansions for the respective critical exponents.) The 
critical exponents of the two systems are the same and 
are given in Table 1 jl^] . Fourier transforming the fields 
V and V in both space and time, the quadratic part of the 
effective action is given by 



oG 

J, 



eff 



?}(-q, -w) [7777 ^(g) - ui4(w)]u(q,w) 



q,a; 

+ ^ / ■vi-q,-(^)u2,o{uj)v{q,uj), (13) 

where mi is trivially zero because we expand around 
the stationary point. In the long wavelength and low 
frequency limit, we approximate rj^^{q) ~ + Kq^ and 
ui,i{uj) « uf f* -I- iwa, where uf^ = -qv' {go) / {1 + -q) , 
go — rjvo + F, and uff" and a are given in Appendix B. 

27r^(^) ((i)')„ - 



The correlation function is 7/2,0(1^) 

\T'h 



The bare propagator is given by Gjjy{q,uj) w 

[—iuja + rjKq^ — X~"^]i where x — [uf f' — 1]^"^ is precisely 
the static response to a monotonically increasing exter- 
nal force calculated in MF theory. At the MF critical 
point = and the bare propagator becomes diffu- 
sive, Gvv{q,i^ = 0) ^ 1/<Z^, while the correlation function 
is static, with Gvv{q,ijj) ^ 5{uj)/q'^. This behavior is the 
same as for the RFIM. In contrast, the bare propagator 
for driven elastic CDWs is diffusive also away from the 
critical point. 

To set up an RG calculation, we need to consider terms 
beyond Gaussian in the action. To obtain an approxi- 
mate effective action, we Fourier transform the fields in 
time, retain only the low frequency limit of the response 
functions, and then transform back to time. In the static 
limit, the ui_2 term has bare value 1/1^2(^1, a;2) ~ w, 
with w ~ ?7^i/"(go)/(l +7?), and Mi,3(a;i, a;2, W3) ^ u — 
rfv"'{gQ)/{l 4- rf). Note that w = at the MF critical 
point. We then obtain 

5eff « S^s + J2l^{^ v.[v.it)]' + ^ vmv.it)?}, (14) 

which is identical to the effective action for the soft spin 
RFIM |l^] , with the local velocities here playing the role 



of the local magnetizations in the RFIM. The generic 
behavior at long time and length scales is therefore iden- 
tical to that of the RFIM. To perform the coarse-graining 
transformation, as usual we integrate out modes in the 
wave vector shell A/6 < g < A (fe > 1), rescale coordi- 
nates as X — hx' and t = V't' and choose the rescaling 
of the fields so that the quadratic part of the action is 
unchanged at the critical point {x~^ = 0). This requires 
z = 2 and yields u'^,^ = 6[-(™+")+2l'*/2+2»u„^„. 

If d > 8, all Umn other than those yielding quadratic 
terms in the action renormalize to zero and become ir- 
relevant at large scales, indicating that the MF line for 
the onset of the vertical jump in the v — F curve (where 
= and w 7^ 0) remains critical The MF scal- 
ing exponents are exact in this case. 



Exponents 


Mean-field 


e expansion 


f} 


1/2 


1/2 - e/6 + ©(e^) 




2 


2-e/3 + 0(e2) 


135 


3/2 


3/2 + 0.0833454e^ 



TABLE I: Universal exponents for the critical point discussed 
in the text. The exponents j3 and 5 describe how iJ scales with 
F and r], respectively, is the correlation length exponent. 

For d < 8 we need to distinguish two cases, depending 
on whether w' ^ h'^~'^/'^w is initially finite or zero. A 
bare action with — and 7^ at the onset line 
Fc{ri) corresponds to the hysteretic region -q > rjc. In 
this case the RG analysis carried out by Dahmen and 
Sethna for the RFIM shows that w flows to 00, suggest- 
ing a sharp jump onset in the hysteresis loop for 77 > rjc- 
A bare action with x~^ = and w = corresponds to 
the critical point of interest here. If the bare w is zero, 
we must consider higher order terms. In particular we 
find that u' ~ b^~'^u, while all the higher order response 
functions are irrelevant for d = 6 — e. So dc = 6 is the 
upper critical dimension and for d > 6 the critical ex- 
ponents near (Fc,r]c) are given by MF theory. Below 6 
dimensions, u is relevant. Corrections to the critical ex- 
ponents were computed in 6 — e dimensions in , the 
results are summarized in Table 1. 

Many open questions remain. First, the RG calcula- 
tion near {Fc,r]c) described here apphes for rj < rjc @. 
More work is needed to extend it above 77c, as well as for 
a proper understanding of the noise in this region. Sec- 
ondly, we would like to study the finite range viscoelastic 
model for general value of r. 

We thank Alan Middleton, Jennifer Schwarz for illumi- 
nating discussions and help with one of the figures. We 
also thank Jim Sethna for helpful comments. MCM was 
supported by NSF through grants DMR97-30678 at Syra- 
cuse and PHY99-07949 at ITP, Santa Barbara. KAD 
acknowledges support from the Materials Computation 
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72783, and an A. P. Sloan fellowship. 



IV. APPENDIX A: MEAN-FIELD SCALING 

Here we derive the scaling behavior near the MF crit- 
ical point. It is convenient to rewrite the mean velocity 
as 

l + r, 

by introducing the function ^{g), 



(15) 



2g 



v{g) = dh p{h) 



(16) 



with p{h) the normalized distribution of pinning 
strengths. The MF critical point, rj — rjc, F — Fc and 
V = Vc, is the point where the mean velocity curve has 
infinite slope. The conditions for the critical point easily 
are expressed in terms v^g) by differentiating both sides 
of Eq. (|l^) with respect to g, to obtain 



(17) 



( 



(18) 



KdF^yr, 

where the prime denotes a derivative with respect to g. 

To determine the critical line Fcij]) where the slope 
of the (v, F) curve diverges, we solve Eqs. ([l5|), ( |l7| ) 
and ( |l^) for rj near rjc- By expanding both sides of 
Eq. (Iia) in Sgdr]) = gdrj) - gc and Sr] = t] - rjc, and 

using that v"{gc) ~ 0, we obtain 5gc{'r]) ^ ± 



-2(5 T) 

^2 ^n' 



with v'd = v"'{gc) < 0. We then use ^^^(77) = dgdri) - 
VcSt] — rjcSvc — SrjSvc, expand both sides of Eq. ( p^ ) and 
eliminate Sgdri) between the two equations, to obtain 



-Sr] ± 



3?7c'(l + ^c) 



1/2, 



{Sr]) 



3/2 



(19) 



Note that the critical line Fc{ri) is rotated as compared 
to the critical line in the magnetic system. 

To obtain the scaling function, we now write v{F, rf) = 
v{Fc,r]c) + Sv in Eq. (p^, expand around the critical 
point and eliminate 5g in favor of Sv, Srj and SF. Re- 
taining only terms up to cubic in the deviations from the 
critical point, we obtain a cubic equation for Sv, given by 

-{SvY ^ SrjSv H — Srj - 



3!(l + ?7c) 



1 



Vc 



1 



SF = 



Vc 



(20) 

We define a rotation from {F, rf) to new scaling variables 
[F' , 77), with F' = F + i^^V- terms of the new vari- 
ables Eq. (pH) becomes 



(Sv)^ 



vyj' 



SrjSv -\r 



6(1 + ^c),^, 



SF' = 0, 



(21) 



whose solution is given by Eq. (^, with Q±{x) the small- 
est real root of the cubic equation 



^3 , 6 6(1 + r;,) 



vlK'' 



vt< 



(22) 



The scaling function satisfies G±{0) = constant and 
g±{x » 1) ~ 



V. APPENDIX B: EVALUATION OF THE 
LINEAR RESPONSE 

The linear response function wi 1 is obtained by per- 
turbing the system with a (5-function pulse of strength 
e > at time ti and evaluating the response (u(i))/.c at 
time t > ti, 



(23) 



where u{t) is obtained by solving Eq. (^ for t{t,ti) = 
eS{t — ti) (to simplify the notation, we drop the sub- 
script on the local field hi Eq. (12)). The other local 
response and correlation functions can be evaluated by 
similar methods. 

As shown in Fig. ||, the perturbation has two effects: it 
yields a discontinuous jump of the displacement at time 
ti and it causes a shift in the "jump time" tj where the 
displacement jumps between neighboring wells of the pe- 
riodic pinning potential. Without loss of generality we 
assume tj <ti < tj + T, where T is the period defined 
in Eq. (||) and tj is the jump time of the unperturbed 
solution. 




o perturbed 
o 



u=2 



u=l 



u=0 



t +2T 

time 



FIG. 3. Time evolution of the displacement u{t). Both the unperturh 
shown. The perturbed solution is obtained by applying a 5-function puls( 

For t >ti, the perturbed solution can be written as 
u{t,ti;tj)^eit~h)u^np{t + Stj;tj) , (24) 
where Uunpit;tj) is the unperturbed solution, given by 
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Munp(i; tj) 



E 

n=0 



e{tj + {n + 1)T - t)Q{t -tj- nT) 



-\(t-tj-nT) 



'XT 



(25) 



with A = + 77). The shift in jump time 5tj is de- if tj > ti — T + t^, with = j ln(l + ea), and 6tj — 
termined by requiring u{tj — Stj,ti;tj) = 0, with the otherwise, with a = ri{e^'^ — 1)/(1 + rj). 
result, The local response is obtained by averaging over a uni- 

form distribution of jump times, as well as over the distri- 
bution p(h) of pinning strengths. The average over jump 
times of the perturbed velocity is given by 



Stj = --ln(l 
A 



-X{tj-\-T~ti)\ 



(26) 



*^-^+*^ dh d r 



ti-T 



T dt 



Q{t - ti)Uunpit-tj) 



dti d 



eit-ti)u^npit + Stj-tj) 



(27) 



Carrying out the integrals p8| , differentiating the response with respect to e and evaluating the result at e = 0, we 
obtain the mean local linear response function, xiit — ^i) = ui^i{t,ti), with 



Xi{t)=S{t) 



V 



i + v 



-^^^e(t^nr)e((n + i)r 
1 + v)T 



t)e 



'\(t~nT) 



n=0 



{nT-t+±;)ie'^-l)-l\)., (28) 



where the brackets (...) denote the average over p{h). For 
t > 0, the response consists of the sum of a (5-function 
contribution and a periodic function. The frequency de- 
pendent local response is given by 



(29) 



Jo 

l+n \T ^ 



where the limit C — > is intended, Us = 2t: s/T, with s 
an integer, and 



ryAe^^ (1 - e-^^y im, 
1 + 1] AT {iujs - A)2 



(30) 



For small frequency, we can carry out the summation 
over s and find 



with 



and 



Xi{uj) « uf f +iuja 



Stat 
"1,1 



l+f] 



(31) 



(32) 



l + r/\A 2A2r 
-^(e--l)(l-e--)) 



(33) 
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